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, The goal to describe plutonium phases from "first principles" calculation methods is complicated 

by the problem of 5f-electrons localization. While for early actinides (Th, U,Np) standard DFT 
(Density Functional Theory) description with itinerant 5f-electrons works well for late actinides 
(Am, Cm) DFT calculations with completely localized (pseudocore) 5f-electrons give satisfactory 
results. However plutonium presents a border case of partial localization and both limits (itinerant 
and completely localized) are not valid. We present a review of the methods used to solve this 
' problem and discuss what could be the reasons for their successes and failures. 
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PACS numbers: 71.27.+a, 71.70.-d, 71.20.-b 



I. INTRODUCTION 



The Density Functional Theory (DFT)i in Local Density Approximation (LDA) (or in more elaborated Generalized 
Gradient Approximation (GGA)) has been very successful in describing equilibrium volumes and stable crystal struc- 
• tures for pure elements and its compounds. While there are sometimes difficulties in reproducing spectral properties, 
such as energy gap values for semiconductors and insulators, which are defined by excited electronic states, ground 
Q state energy as a function of the crystal volume and atomic positions is reproduced usually quite well. For rare earth 
elements problems appears in DFT calculations when one treats 4f-electrons on equal footing with all others. However 
in this case a good approximation is to consider 4f states as completely localized in the same way as core orbitals (so 
I 1 called pseudocore approximation). 
, Actinide elements give a challenge to DFT because 5f-electrons are not obviously localized as 4f electrons in rare 
Q\ ' earth but also can not be always considered as itinerant. On the Fig^ experimental values of equilibrium volume for 
5d, 4f and 5f elements are presented. One can see that 5d elements show decreasing of the volume as number of 5d 
electrons increases till 5d shell becomes half filled. Further filling of 5d shell results in increasing of the volume. That 
corresponds to the strong participation of 5d electrons in the chemical bonding of transition metals evidencing their 
fully itinerant behavior suitable for treatment by DFT. In contrast to that lanthanides show very weak dependence 
of the volume on 4f shell filling (with the exception of Eu and Yb case where valency state is +2 instead of the 
common for lanthanides +3). That fact can be explained by fully localized nature of 4f electrons and hence absence 
of their contribution to chemical bonding. In actinides series beginning of the curve follows transition metals pattern: 
decreasing of the volume with increased filling of 5f-shell from Th till Np. However with further increasing of the 
I ■ number of 5f-electrons there is a sizable jump of the volume value for delta-Pu and Am and for late actinides Cm,Bk 
and Cf one observes weak dependence of the volume on 5f-shcll filling resembling lanthanides with localized 4f electrons. 
This curve can be understood if one assumes itinerant nature of 5f-electrons from Th to Np, localized 5f states from 
Am to Cf and intermediate situation that can be described as "partial 5f-electrons localization" for PvA The problem 
is complicated by the fact that Pu itself exists in 6 various allotropies with volume values differences reaching 20% 
(e.g. between alpha and delta phases). That can be interpreted as more localized 5f-electrons in delta phase than in 
alpha phase. 



II. STANDARD DENSITY FUNCTIONAL CALCULATIONS 



Standard DFT calculations for experimentally observed paramagnetic ground state gave good results for equilibrium 
volumes of early actinides from Th till Np (Fig[2Jl. However already for alpha phase there is an underestimation of 
theoretical volume comparing with experimental one and for delta phase the disagreement is rather large (more than 
20 %). These results agree well with the above analysis of the experimental data for equilibrium volume of actinides. In 
DFT all electrons are considered to be itinerant and participating in chemical bonding. That is a good approximation 
for early actinides but an obviuosly wrong one for late actinides starting from Am and Pu in delta phase, where 
5f-electrons are completely or partially localized. 
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The term "partially localized" means that electrons on every one of 5f-orbitals can demonstrate properties char- 
acteristic for localized and itinerant pattern. That can be understood as in different periods of time electrons can 
either be sitting on specific atomic sites or being spread over the crystal. This effect can be described by Dynam- 
ical Mean-Field Theory (DMFT)2 with time or energy dependent self-energy operator. In so called "mixed levels" 
scheme of O. Eriksson et al» "partial localization" was imitated by treating some specific 5f-orbitals fully localized 
(pseudocore) and all others fully itinerant. It was found that the best agreement for equilibrium volume value for Pu 
in delta phase was obtained (Fig|3J) with four 5f-electrons considered localized. Similar approach was used in so called 
Self Interaction Correction (SIC) calculation scheme where calculations gave the best description of Pu in delta phase 
with localization of three 5f-electrons£. 

Recently it was found that combination of Generalized Gradient Approximation (GGA)) and spin-polarization taken 
into account results in drastic improvement of DFT results for equilibrium volume values§*£*£. While for alpha phase 
non-magnetic and magnetic solutions gave total energy as a function of volume curves with close minimum positions, 
the delta phase calculations show sizable increasing of calculated equilibrium volume value when spin polarization was 
taken into account (Fig0J. Not only delta phase but also other crystal structure phase have their calculated volume 
values in much better agreement with experiment (FiglSJ) . The reason for such effect can be understood by observing 
that phases with large crystal volume have much larger magnetic moments values (Fig|SJ). While for low volume 
alpha phase spin magnetic moment values per Pu ion are « 2/ljb, in large volume delta phase the corresponding value 
is ~ 4:[Ib 6 - Spin polarization results in exchange splitting of 5f-orbitals and hence partial removal of 5f states from 
the Fermi level. That leads to suppression of 5f contribution to the chemical bonding and increasing of calculated 
equilibrium volume values. One can say that strong spin polarization leads to "partial localization" of 5f-clectrons. 
This effect is most pronounced in the case of Curium^ where configuration f 7 results in full spin polarization with 
S = 7/2 and maximum value of exchange splitting for 5f-orbitals with complete removal of 5f states from the Fermi 
energy. 

The price for this improvement of DFT results is magnetic ground state of Plutonium obtained in calculations in 
disagreement with experiment. In the recent review of J.C.Lashley et alifi it was shown that not only there is no 
experimental evidence for long range magnetic ordering but also the possibility of disordered local moments presence in 
all Pu phases can be excluded by the results of whole set experimental data including magnetic susceptibility, NMR, 
specific heat in magnetic field, neutron elastic and inelastic scattering. Recent results^ of muon spin relaxation 
measurements on elemental Pu set an upper limit on ordered moments for alpha-Pu and Ga-stabilized delta-Pu at 
T=4K of 0.001 fiB- This disagreement of spin-polarized DFT calculation results with magnetic measurements can not 
be explained by cancellation of spin magnetic moments by the antiparallel orbital moment 12 . The spatial dependence 
of the orbital and spin magnetizations is different around the nucleus and if the total moment is equal to zero the 
difference in their spatial extent would still allow a measurable signal to be seen in neutron scatteringiS. Also the 
arguments for magnetic moment cancellation ini£ were based on using formula A/ iot = hb{L — 2S) that gives zero 
with L = 5 and S = 5/2. However formal total moment value for such case is J = 5/2. 

In order to understand the origin of spin polarized DFT calculation problems it is instructive to consider much 
simpler case than Pu with its "partial localization" and complicated phase diagram. Americium at ambient pressure 
has simple hexagonal close-packed crystal structure with localized 5f-electrons in / 6 configuration in jj coupling 
scheme with zero values of spin, orbital and total moment. However DFT calculations^ gave fully spin polarized 
solution with huge values of spin and total moments (Fig.JTJHJ))- It is not only magnetic properties that DFT gave 
wrong but also there is a strong disagreement with spectroscopy experiments. In Am photoemission spectra (Fig|5J) 
occupied 5f band is centered around 3 eV, while in calculated density of states (Fig.©) its position is only 1 eV 
below the Fermi energy. This fact is consequence of full spin-polarization when 5f spin- up sub-band (Fig. d>) with 
full capacity 7 electrons is filled with only 6 and Fermi level is inside spin-up sub-band. The ground state obtained 
in DFT calculations corresponds to the LS coupling scheme with maximum possible spin moment value. In contrast 
to that in jj coupling scheme the good quantum number is only total moment J and there are j = 5/2 sub-shell 
with capacity 6 electrons and j = 7/2 sub-shell with 8. Then Am with / 6 configuration has fully occupied j = 5/2 
sub-shell and empty j — 7/2 sub-shell with J = 0. 

The problem of which coupling scheme (LS, jj or intermediate) will be realized for particular ion is defined by the 
competition between Coulomb exchange interaction preferring maximum spin polarization with strong energy splitting 
between spin- up and spin-down subbands and spin-orbit coupling leading to energy separation between j = 5/2 and 
j = 7/2 subsells. While for 3d ions spin-orbit coupling is weak and LS coupling scheme is always valid, for 5f-orbitals 
of actinides spin-orbit coupling is strong enough to compete with exchange interaction. In the result for Am with / 6 
configuration jj coupling is realized with S = L = J = and for Cm with f 7 configuration LS coupling wins with 
S = 7/2. The failure of DFT in describing Am can be due to the overestimation of exchange interaction strength. In 
DFT the functional depends on spin-up and spin-down electron densities. This form of the functional suggests that 
spin is a good quantum number while for the case of strong spin orbit coupling it can be not true any more. 
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III. COULOMB INTERACTION EFFECTS IN STATIC MEAN-FIELD APPROXIMATION: LDA+U 

CALCULATIONS 

The problem of "first principles" calculations of electronic structure and ground state properties of Plutonium 
is determined by the question how to describe 5f-electrons localization. The physical origin for localization are 
correlation effects due to Coulomb interaction between 5f-elcctrons. In order to include correlation effects into DFT 
calculations new methods were developed. Two of them are LDA+Ui4 and LDA+DMFTiSJi methods. In the first 
one Coulomb interaction is treated in static mean-filed approximation (unrestricted Hartree-Fock) and the second one 
via Dynamical Mean-Field Theory (DMFT) with energy (or time) dependent complex self-energy operator. 

In LDA+U method potential is orbital dependent with occupied orbitals having lower energy than unoccupied ones. 
In the result 5f-orbitals are moved away from the Fermi energy and their contribution to chemical bonding is strongly 
reduced comparing with paramagnetic DFT results. LDA+U calculations for Pu in delta phas e - 12 ! 15 gave significant 
increasing of equilibrium volume value to a good agreement with experiment (Fig.^UJ). However in these works the 
same problem as in spin-polarized DFT calculations appears: strong spin-polarization with a large values of magnetic 
moments. 

Recently it was found that LDA+U equations can give a nonmagnetic solution 18,19 for Pu with 5f-shell in ground 
state with S = L = J = and calculated equilibrium volume for delta Pu in a good agreement with experimental 
value. The physical origin for such ground state is strong spin-orbit splitting of 5f states that is larger than 5f band 
width. In pure LDA calculations without LDA+U potential correction (Fig.l|llh,)) 5f-orbitals density of states (DOS) 
consists of two well separated sub-bands being formed predominantly by orbitals with total moment j = 5/2 for the 
occupied states and j = 7/2 for the empty ones. Such 5f-DOS pattern means that nonmagnetic state with filled 
j = 5/2 sub-shell is already nearly preformed in standard LDA description. There is however strong admixture of 
j = 5/2 states to the empty band and j — 7/2 states to the occupied one. When LDA+U potential correction is 
applied occupied orbitals are shifted down in energy and empty ones up. In the result fFig. l|llb )) one obtains pure 
j = 5/2 and j = 7/2 sub-bands with increased energy separation between them. The Fermi energy is positioned 
on the top of the occupied j — 5/2 sub-band that agrees well with the peak position around 1 eV in photoemission 
spectra (Fig. Q )- 

LDA+U calculations for Am22*2I gave a nonmagnetic solution with 5f-shell in ground state with S = L = J = 0. 
The LDA+U results for Am are in much better agreement with experiment than spin polarized DFT calculations 
not only in magnetic but also in spectral properties. The general 5f-DOS pattern fFig. l|13ll4|l % ) is the same as for 
Pu fFig. (|ll|l ) with the important difference that position of the Fermi level is not on the top of occupied j = 5/2 
sub-band as it was the case for Pu but at the bottom of empty j = 7/2 sub-band and smaller 5f band width due 
to larger volume of Am. This shift of the Fermi level is due to additional valence electron in Am comparing with 
Pu which goes into s— ,p — d— states not shown on Fig. (|I3ll4|) . In calculated 5f-DOS occupied 5f band is centered 
around 3 eV in a good agreement with Am photoemission spectra (Fig|5J|. 

While LDA+U method solved the problem of 5f-electrons localization without developing magnetism that was not 
observed in experiment there are still important disagreements of LDA+U results with experimental data. First of 
all it is the absence in the LDA+U calculated spectra of the sharp peak on the Fermi level observed in experimental 
photoemission spectra fFig. lfL^)) ). While peak at 1 eV in calculated spectrum for delta Pu agrees well with the 
corresponding feature in the experimental spectrum, the sharp peak close to the Fermi energy in photoemission 
spectrum does not find any correspondence in calculations. 

IV. COULOMB INTERACTION EFFECTS IN DYNAMIC MEAN-FIELD THEORY: LDA+DMFT 

CALCULATIONS 

This fact can be understood as a manifestation of the "partially localized" nature of 5f-electrons in Plutonium. 
In fully localized case one has occupied lower Hubbard band below the Fermi energy and empty upper Hubbard 
band above it. If electrons are itinerant then partially filled quasiparticle band crossing the Fermi energy is a correct 
description. For the intermediate case of "partial localization" both features: lower Hubbard band below the Fermi 
energy and partially filled quasiparticle band would be present in the spectral function, with additional effect of 
quasiparticle band narrowing due to correlation effects. That effect can be reproduced in Dynamical Mean-Field 
Theory (DMFT) 3 with its famous "three feature structure" in the spectral function (lower and upper Hubbard bands 
with quasiparticle peak on the Fermi energy between them). LDA+U method as static mean- field theory can give a 
good approximation for the Hubbard bands but not the quasiparticle peak. This peak needs dynamical fluctuations 
to be taken into account for proper description. 

Another deficiency of LDA+U solution for Pu is 5/ 6 configuration of the calculated ground state that corresponds 
to completely filled j = 5/2 sub-band. From analysis of absorption spectra of Pu2£ mainly 5/ 5 configuration can be 
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found. The reason for this disagreement could be the following. The Fermi level in LDA+U solution for delta Pu 
(Fig.ljnj))) is positioned on the top of occupied j — 5/2 band with the center of the band less than 1 eV below the 
Fermi energy. That can result in high probability of 5f-electrons excitation from occupied j = 5/2 band to the empty 
s— ,p — d— states above the Fermi energy. This is dynamical fluctuations effect and if it will be taken into account 
the average number of 5f-electrons can be decreased from the static mean-filed solution value. 

The disagreement with experiment of LDA+U solution discussed above leads us to the conclusion that in order do 
have correct description of 5f states of Pu one should explicitly take into account dynamical fluctuations in calculations. 
That can be done only in LDA+DMFT methocUSiii. The essence of DMFT is mapping of the problem for lattice 
of atoms with Coulomb interacted electrons on the effective impurity problem for an ion interacting with effective 
bath (reservoir) characterized by energy (time) dependent hybridization function that is calculated self-consistently 3 . 
This effective impurity problem should be solved taking into account full Coulomb correlations between electrons on 
the ion in reservoir. While being simpler than full lattice problem the impurity problem can be still very expensive 
computationally, especially for multi orbital case. 

There were developed many methods to solve effective impurity problem, so called "impurity solvers" (for the 
recent review on this problem see^). One group of them form approximated methods either based on perturbation 
theories in hybridization strength as Noncrossing Approximation (NCA) or in Coulomb interaction as FLEX. To the 
same group belongs interpolative approaches like IPT and its extensions. Another group consists of the methods 
that can be considered formally exact: Quantum Monte Carlo (QMC), Exact Diagonalization (ED) and Numerical 
Renormalization Group (NRG). The problem with approximated methods is that they give uncontrollable accuracy 
results. Even if those methods were tested by comparison of their results with exact approaches for some simple cases 
it is no guaranty that they will work for the problem with a large orbital degeneracy. ED and NRG methods are not 
practical for orbital degeneracy larger than two. The only method that can in principle give reliable and practically 
achievable solution is QMC. However for 5f-orbitals with degeneracy equal seven full DMFT-QMC solution can require 
use of the most powerful computers. 

LDA+DMFT methodi£*i£ was applied to Pu problem^ and encouraging results were obtained including the possi- 
bility of the double minima in energy vs volume curve (Tig- Hf 5(1 ) and the peak on the Fermi level fFig. (|I 6|) L However 
for "impurity solver" the authors used interpolative approach with a simple analytical form for self-energy with 
parameters being adjusted to obey known various asymptotes. 

Another attempt to apply LDA+DMFT method for Pu was done in2£ where the authors had started from nonmag- 
netic LDA+U solution and included fluctuation via "Spin-orbit T-matrix FLEX approach" based on the perturbation 
theory in Coulomb interaction strength parameter U. The calculated spectral function (Fig. I^TJ) ) shows in comparison 
with experimental spectrum too strong intensity of the peak on the Fermi energy and suppressed lower Hubbard band. 

Recently we have done LDA+DMFT calculations for Pu in delta phase with QMC "impurity solver"— where 
only j = 5/2 5f-orbitals were treated as fully dynamical while j = 7/2 5f-orbitals were described by static mean- 
filed (Hartree-Fock) approximation. The justification for this was LDA+U solution (Fig. (|f f b )) where unoccupied 
j = 7/2 band was found to be situated well above the Fermi energy. The calculated spectral function (Fig.^HJ)) 
shows correct position of the both features of the experimental photoemission spectra: lower Hubbard band at « 1 
eV and quasiparticle peak on the Fermi level but the relative intensity of quasiparticle peak was underestimated. The 
average number of 5f-electrons due to the dynamical fluctuations was significantly decreased from its static mien-filed 
(LDA+U) value of six electrons to «5.5 that corresponds to equal weights of 5/ 5 and 5/ 6 configurations in the ground 
state. 



V. CONCLUSION 

In conclusion we have shown that the source of the problem in describing Pu from "ab-initio" electronic structure 
calculations is "partial localization" of 5f-electrons. Both limits, itinerant as in DFT and completely localized treating 
5f-electrons as pseudocore are not appropriate for Plutonium problem. "Partial localization" means that 5f-electrons 
part of the time spend sitting on particular Pu ion and the rest of the time are spread over the crystal. That 
corresponds to dynamical (time dependent) fluctuations which can be described by Dynamical Mean-Field Theory 
(DMFT). Effective impurity problem appearing in DMFT must be solved by Quantum Monet Carlo (QMC) as the 
only reliable and practically realizable method. However that will require heavy use of modern multi processor 
computers. 
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FIG. 1: Experimental values of equilibrium volume for 5d, 4f and 5f elements as a function of number of electrons in partially 
filled shelM. 
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FIG. 2: Theoretical values of equilibrium volume for 5f elements from Density Functional Theory calculations 27 . 
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FIG. 3: Total energy as a function of volume for Pu in delta phase from DFT calculations with various number of 5f-electrons 
treated as completely localized ("mixed level" scheme O. Eriksson et al(-) 
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FIG. 4: Total energy as a function of volume for alpha and delta plutonium from spin-polarized DFT calculations of A.Kutepov 
et al£. 
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FIG. 5: Total energy as a function of volume for various crystal structure phases of plutonium from spin-polarized DFT 
calculations of P.Soderlind et al-. 
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FIG. 6: Magnetic moment for Pu in various crystal structure phases from spin-polarized LDA calculations (P. Soderlind et al-) 




FIG. 7: 5f-electrons partial density of states for Am from spin-polarized LDA calculations (P. Soderlind 
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FIG. 8: Magnetic moment for Am with spin and orbital contributions from spin-polarized LDA calculations (P. Soderlind et 
all!) 
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FIG. 9: Experimental photoemission spectra of Am, Pu and Sm (J.R.Naegele et a^) 
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FIG. 10: Total energy as a function of volume for delta phase of Pu from DFT and LDA+U calculations ( J.Bouchet et al— ) 
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FIG. 11: Partial 5f densities of states for Pu in delta phase from LDA (a) and LDA+U (b) calculations (A.Shorikov et al— ) 
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FIG. 12: Comparison of experimental and calculated (LDA+U) spectra for alpha and delta Pu and PuSb (A.Shorikov et al— ^) 
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FIG. 13: 5f-electrons partial density of states for Am from LDA calculations (A.Shorikov et 
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FIG. 14: 5f-electrons partial density of states for Am from LDA+U calculations (A.Shorikov et 
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FIG. 17: Comparison of experimental and calculated (LDA+U, DMFT-FLEX) photoemission spectra for delta Pu (L.V. 
Pourovskii et al— ). 
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FIG. 18: Comparison of experimental and calculated (DMFT-QMC) photoemission spectra for delta Pu (J.Kunes et al— ) 



